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We show that the competition between interactions on different length scales, as relevant for 
the formation of stripes in doped Mott insulators, can cause a glass transition in a system with 
no explicitly quenched disorder. We analytically determine a universal criterion for the emergence 
of an exponentially large number of metastable configurations that leads to a finite configurational 
entropy and a landscape dominated viscous flow. We demonstrate that glassines is unambiguously 
tied to a new length scale which characterizes the typical length over which defects and imperfections 
in the stripe pattern are allowed to wander over long times. 
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I. INTRODUCTION 

The competition of interactions on different length 
scales is one of the mechanisms able to stabilize mesoscale 
phase separations and create spatial inhomogeneities in 
a wide variety of systems. The most typical situation 
is a competition between short-ranged forces, that fa- 
vors the formation of a uniform condensed phase, and 
the long-range forces which can energetically frustrate 
this condensation. Classical examples are the formation 
of domains in magnetic multilayer compoundsEra, meso- 
scopic structures built by assembling nnalymers in solution 
or amphiphilcs in water-oil mixturesau.Very often, these 
systems exhibit a long time dynamics similar to the re- 
laxation seen in glasses. Conversely, many proposals have 

been made that the glassy behavior of molecular liquids 
might arise from frustration of specific crystalline- like or- 
ders incompatible with global packing, e.g., icosahedral 
order in dense liquids^ . 

The observation of complex orbital and charge pat- 
terns in CMR-mangaiiites or of charge stripes in doped 
nickelates and cupratea3'Ll suggests that a similar compe- 
tition causes inhomogeneoua structures in these strongly 
correlated electron systemsB. This point of view is Haft, 
ported by the observation of spatial inhomogeneitiesETcJ 
as well as slow, activated glassy dynamicsEM, as seen in 
recent NMR experiments. Upon cooling, the Cu-NMR 
signal in La2- x Sr x CuOi based systems disappears (is 
"wiped out"). This has been interpreted in terms of an 
electronic relaxation-slower than the Larmor precession 
of the nuclear spinsHL-l As a result, the NMR signal 
decays so fast that it simply cannot be detected an** 
more. Thus the "wipe out effect" discussed in Refs.EliL3 
is clear evidence for a dramatic increase of the relax- 
ation times of the electronic system. Similarly, La- NMR 
was used to directly show that there is a glassy acti- 
vated dynamics with a maximum of T^ 1 (T), separating 
relaxational dynamics which is slower than the nuclear 



Lamor frequency . at l ow temperatures from faster pro- 
cesses at higher Tt3ta. The typical activation energies of 
these|-oynamical processes have been analyzed by Curro 
et al.li-3 who made the surprising observation that they 
are rather independent of the specific details of added 
impurities etc. Also, the width of the distribution of 
activation energies is comparable to its mean value in 
systems with rather different chemical composition. In 
view of this striking universality of the anomalous long 
time relaxation in doped Mott insulators, we recently 
suggested that glassiness in these systems is self gener- 
ated, i.ep-,it| does not rely on the presence of quenched 
disordenlSO. The latter may, in general, further sta- 
bilize a glassy state. This is supported by molecular 
dynamics calculations for charge ordering in transition 
metal oxides, which found an anomalous long timepcelax- 
ation with a power spectrum similai_.to 1 / /-noisecS. In 
addition, cecently M^rkiewicz et alS3 analyzed neutron 
diffractionEj, NQRO, ^SREL anelastic relaxationEil and 
susceptibility measurementsEj, spanning altogether more 
than 10 orders of magnitude of frequency, and also find a 
"universal" behavior of the activation energies in under- 
doped cuprates. They also find a good description of the 
relaxational dynamics using a Vogel-Fulcher law which 
we predicted based on a entropic droplet argumentEj. 



In Ref.ESEZI we showed that the competition between 
interactions on different length scales causes the emer- 
gence of an exponentially large number of metastable 
states (with the system size). This is generally consid- 
ered as a condition for the anomalous dynamical features 
of glassiness, like aging, memory effects and ergodicity 
breaking. It is also,the heart of the random first order 
transition scenarios for vitrification of molecular liquids, 
originally motivated by the similarities Jaetween density 
functional theories of aperiodic crystalsc3 and the mean 
field theories for random spin-glasses. This scenario is 
now believed to apply to a much more general class of 
systems. Our result was obtained using a new replica 
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approach, EaEZl and by solving the resulting many body 
problem numerically within the self consistent screening 
approximation. In this paper we develop an analytical 
approach which enables us to identify the underlying 
physical mechanism for glassiness in a uniformly frus- 
trated system. We furthermore discuss that our results 
can also be obtained within a dynamical approach, where 
glassiness is associated with an unconventional long time 
limit of the charge correlation function. 

In the next section we introduce the model we investi- 
gate and summarize the main results of this paper. The 
details of our approach are presented in section IV, sub- 
sequent to our summary of the aspects of the stripe liquid 
state that will be important for our results in section III. 
In section V we conclude and give a list of further open 
questions. 



II. MODEL AND OVERVIEW 



This paper develops an analytical approach to glassi- 
ness in a uniformly frustrated system. By uniformly frus- 
trated we mean that there is a competition of interac- 
tions on different length scales and there are no explic- 
itly quenched degrees of freedom, like the ones caused by 
additional defects or imperfections. We study a model 
with local tendency towards phase separation, frustrated 
by a long range interaction, which, as we will show, has 
all necessary features to exhibit a glass transition and yet 
is simple enough to be treated analytically. In the con- 
text of cuprate systems the model has been proposed by 
Emery and Kivelsontl and is defined by the Hamiltonian: 

7i= l -jd d x {r ^(x) 2 + (V^(x)) 2 + ^(x) 4 } 

Here, y(x) characterizes charge degrees of freedom, with 
</?(x) > in a hole-rich region, </?(x) < in a hole poor 
region, and <p{x) = if the local density equals the av- 
eraged one. If r < the system tends to phase sepa- 
rate since we have to guarantee charge neutrality^) = 0. 
The coupling constant, Q, is a measure for the frustra- 
tion between this short range coupling and the long range 
Coulomb interaction. For Q = and ro < we expect at 
low temperatures long range charge ordering. The or- 
dering temperature can be estimated within mean field 
theory as T° = 2,r , with momentum cut off A of the 
order of an inverse lattice constant. As shown in Ref.Ea, 
within a large N approach (where y(x) is generalized 
to an V-component field), for all Q > , the Coulomb 
interaction suppresses this ordered state at finite T. In- 
stead, as revealed by a mean field analysis of Eq.[l] where 
N — > co, the system undergoes several crossovers. 
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FIG. 1. The competition between different lenght scales 

according to the mean field sollution of Eq.[l] 

As can be seen in Fig.[j], at high T two characteristic 
length scales occur. One is the charge ordering correla- 
tion length, £, and the other, Id ~ £ _1 Q -1 ' 2 , is the ef- 
fective Debye screening length of charged regions of size 
£ (see appendix C) . Here we can already recognize the ef- 
fect of the competition between the short-range ordering 
interaction and the long-range Coulomb interaction: the 
charge density is homogeneous within regions of size £, 
but behaves like a plasma, with screening length Id ^> £ 
on larger scales. In Ref.ca it was argued that the emer- 
gence of the screening length Id supports the formation 
of compact ordered domains of size Id which then give 
rise to a slow motion and thus glassiness. How such De- 
bye screening should cause such compact domains and 
glassiness has not been made explicit, however. Our 
replica approach gives no indication for glassiness in the 
temperature regime where the Debye screening theory 
applies. Thus, despite using the same Hamiltonian, the 
stripe glass phase discussed in this paper is qualitatively 
different from the scenario of RefO. A more detailed 
discussion of these aspects is given in appendix C. 

As the temperature is lowered, £ increases monotoni- 
cally while Id decreases until it becomes of the order of £. 
At this point, the Debye screening approximation breaks 
down and the system crosses over to a new regime char- 
acterized by spatial charge modulations, called stripes, 
with period l m ~ 2irQ~ 1 / 4 and coherence length £tB. 
These modulations are particularly relevant at low tem- 
peratures where the correlation length £ is larger than 
the inter-stripe distance l m and, as we will show, where 
the stripe glass phase emerges. 

The charge correlation function Q (x,x') = 
T _1 (ip (x) ip (x')) of the liquid state at low temperatures 
is then characterized by two length scales, l m and £. Af- 
ter transformation into momentum space, this function 
obeys the following scaling behavior (neglecting effects 
due to anomalous powers): 

G(<i)=Cg(ql m ,lm/0, (2) 

with Q (q) peaked at the modulation wave vector q m = 
with broadening In thermodynamic equilib- 

rium, the model, Eq.[l] undergoes a stripe liquid - stripe 
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solid transition at same temperature T c . Within a_spher- 
ical approximationEj or a large- iV approximationcj, the 
transition is of second order and T c — > (unless one takes 
additional lattice corrections into account which yield a 
finite Tjia). In the Ising limit (N = 1), T c > and there 
are indications that the transition is driven first order by 
fluctuationscil. 

Our results indicate that this phase transition may not 
be reached kinematically in which case the system under- 
goes a glass transition instead. In fact, according to our 
theory, glassiness emerges if the inter-stripe correlations 
in the stripe liquid phase are sufficiently strong; specifi- 
cally if the ratio £/l m is larger than a critical value which 
we find to be close to 2. The temperature, Ta, where this 
happens, within large- N approximation of Eq.jj], is given 
by 



Due to this additional length scale the following scaling 
behavior of the long time correlations results: 



(3) 



which decreases for increasing frustration parameter Q . 

Note that the criterion f/Z m ~ 2 for glassiness is likely 
to be much more general than the specific formula for Ta, 
which depends on details of the model. Also, in a more 
realistic model, the stripe liquid - stripe solid transition 
is expected to be of first order due to an additional term 
~ ip 3 in Eq.|l| which exists if particle hole symmetry is 
broken. Nevertheless, our results do not depend on an 
actual divergence of £ but solely that it is larger than a 
few inter-stripe separations. Therefore, our theory also 
applies if the transition is only moderately first order and 
the stripe solid does not occurs unless £ > 2l m . In this 
case, the equilibrium transition is avoided and a glassy 
state results. 

Our theory yields that below Ta the system estab- 
lishes an exponentially large number of metastable states 
and long time correlations, characterized by the correla- 
tion function !F (x, x') = T _1 lim t _, 00 (</> (x, t) <p (x', 0)). 
These long time correlations occur even though no state 
with actual long range order exists. Even more interest- 
ingly, long time correlations with T (x, x') ^ are unam- 
biguously tied to a new length scale, A, which character- 
izes the typical length over which defects and imperfec- 
tions in the stripe pattern are allowed to wander over long 
times. This length can be associated with the allowed vi- 
brational motions in a potential minima of the complex 
energy landscape of the system. In analogy with struc- 
tural glasses we therefore call it the Lindemann length of 
the stripe glass. 

Evidently, in the liquid state A is infinite. In the glassy 
state we find that A jumps discontinuously to a finite 
value Xa — £ (Ta) /3 and continuously decreases at lower 
temperature. The discontinuous jump in A, unaccompa- 
nied by a latent heat, indicates the transition is of the 
random first order typeE-3, i.e., has one step replica sym- 
metry breaking. Thus, the glassy state is stable only if 
the slow motion of glassy textures is confined to a range 
smaller than A^ (justifying our term Lindemann- length). 



■ T I • I- ' / I ^ m ^™ 



(4) 



with j (x, y,z) = g (x, y) - g (x, ^ z 2 ~ y 2 ) and g from 
Eq|. 
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FIG. 2. The Q and T correlation functions as given by 
Eqs. | and § for l -f = tt/10 and if- = 4n/5. 

In Fig.^| we show the momentum dependence of T (q) 
in comparison with £7(q). Close to q m = 2ir/l m (if 
||q| -3ml < £ _1 ) w e have F (q) ~ G{q)- Configura- 
tions which are close to perfect stripe arrangements are 
solely characterized by a momentum independent Debye- 
Waller factor, such that T (q) ~ ^ ( c l)- ( - )n tn c 

other hand, if ||q| — q m \ > A -1 long time correlations 
are much reduced compared to instantaneous correlations 
with T (q) ~ A~ 2 (? (q) 2 - These 'tails' of the correlation 
functions are obviously built up by configurations with 
defects and imperfections of the perfect stripe arrange- 
ment. Therefore, A has to be interpreted as the length 
scale over which defects of the stripe pattern are allowed 
to wander after a long time. The glassy state can only 
be supported if A < \a, it melts if defects are allowed to 
wander too far. 

Glassiness, including a viscous, energy-landscape dom- 
inated long time relaxation, sets in due to the occur- 
rence of pexponentially many metastable states M ms oc 
exp (S c ) -l3 We find that the configurational entropy, 



S c {Ta) ~ Q 3/4 V ~ l m 3 V, 



(5) 



where V is the volume of the system. The shorter the 
modulation length, the larger is the number of possible 
metastable states, which is plausible for simple geomet- 
rical reasons. This clearly demonstrates that locally the 
stripe correlations stay intact in all these configurations. 
Furthermore, it is the packing of stripes with different 
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orientation and the arrangement of defects that distin- 
guishes the many different metastable states. 

In the laboratory, the system will freeze into a glass, at 
some temperature Tq < Ta which depends on the cool- 
ing rate. While this glass transition is purely dynamical 
and distinct from a conventional phase transition, a-Jcey 
feature of the ideal glass transition scenario of Ref.c-3 is 
that the slowing arises from proximity to an underlying 
random first order transition at Tk < Tq, where the con- 
figurational entropy vanishes like S C (T) oc T — Tk ■ Our 
theory gives exactly this behavior with 

S c (T)=I m 3 *(^,I m /() V. 

We find W (s, t) = (1 + s _1 ) (1 - tf + ln(l - (1 - tf) 
which vanishes linearly at a temperature Tk- Below Tk 
the system freezes into an amorphous solid state due to 
this 'entropy crisisEj, even for an infinitely slow cooling 
rate. Usually, the reason for the system to prefer the liq- 
uid state over the solid is entropic. If S c — > there is no 
entropic advantage anymore to be in the liquid state and 
the amorphous solid results, even in equilibrium. 

Freezing into a glassy state below Ta implies that 
within our replica approach the barriers between different 
metastable states are infinite. This however is a conse- 
quence of the mean field character of the replica tech- 
nique. Following Refo we argue that the formation of 
a mosaic pattern with dynamically defined droplets of 
size R of different metastable states will occur. Entropy 
driven transitions between different states lead to dy- 
namical processes with relaxation time obeying a Vogel- 
Fulcher law: 
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(6) 



(7) 



determined by the configurational entropy as well as the 
bare surface tension of entropic droplets cto (<?)■ Finally, 
a simple estimate for <jq, based on a variational argu- 
ment, gives D (Q) oc y/Q to a good approximation. This 
was recently found in numerical simulations of the lattice 
version of Eq.|l| by Grousson et|-a|.E3 

Recently, Markiewicz et alxQ analyzed various ex- 
periments performed on La2- x Sr x CuOi, spanning al- 
together 13 orders of magnitude of frequency, and also 
found a "universal" behavior of the activation energies 
in underdoped cuprates. Interestingly, a good descrip- 
tion of the relaxational dynamics using a Vogel-Fulchcr 
law, Eq.|[ is possible. The analysis of Ref.Ej also yields 
that Tk decreases with increasing doping concentration 
x. Using the relation, l m ~ a (2x)~ J23 between inter- 
stripe distance, l m , and the doping concentration, as well 
as lm ~ 27rQ~ 1 / 4 enables us to determine the doping de- 
pendence of Tk and show that it is properly described 



within our theory. This is shown in Fig.|| in comparison 
with the results as deduced from experiment in Ref.E3. 
Indeed, our theory gives the proper doping dependence 
of Tk- If we neglect the differences between Tk and Ta, 
(see inset of Fig.H below), an approximate formula for 

the doping dependence of Tk is: Tk — 47T a x+1 ■ The typ- 
ical doping concentration on which changes in the glass 
transition temperature occur is xq — =L^ ~ 0.008. Note 
that the experiments analyzed in RefJi3 are solely sensi- 
tive to the spin excitations of the system. We argue that 
due to strong coupling between charge and spin degrees 
of freedom (which is evident from the formation of phase 
or anti-phase domain walls), glassiness of the charge den- 
sity causes the observed anomalous long time dynamics 
in the spin channel. 
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FIG. 3. Comparison of Tk with the experimental data 
analysed by Markiewicz et al. 19 



III. THE STRIPE LIQUID 

In this section we summarize the main results for the 
stripe liquid state needed for our subsequent calculations 
in the glass state. We will mostly use results obtained 
within a the leading contribution of a 1/N expansion. 
We calculate higher 1 /N corrections and show that they 
are small in the low temperature regime discussed herej 
Note also that the numerical solution presented in Ref.Ej 
did take higher 1/N corrections systematically into ac- 
count. The agreement between the results of this chapter 
with these numerical results also supports the neglect of 
1/7V corrections for the stripe liquid state. The latter 
will become essential, however, in the glass state. 

The mean field equations of the stripe liquid including 
1/N corrections have been discussed in detail in Ref.ES. 
Some of these results have already been summarized 
above. In what follows we will mostly discuss the low 
temperature regime where the correlation length, £, ex- 
ceeds the modulation length, l m , of the system. This 
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means we will be concentrated in the high temperature 
region of the Figjl} 

In the mean field approximation the correlation func- 
tion is given by: 



r + q 2 + % 



(8) 



where the parameter r = ro + u (4> 2 ) must be deter- 
mined self consistently. At high temperatures, r > 2^[Q 
and the system is characterized by a correlation length 
£ r*j r^ 1 / 2 , similar to the unfrustrated system, as well 
as a Debye screening length Id ~ S.^Q 1 ^ 2 characteriz- 
ing conventional screening of charged objects with linear 
size, £ and charge ~ Q 1 ^ 2 . At lower temperatures, when 
r(T) < 2\/Q, simple Debye screening breaks down and 
the system establishes modulated structures with modu- 
lation length (inter stripe distance), l m — 4ttJ \J2\fQ — 
and correlation length , £ = 2/\/r + 2v^.El 



2/^r + 2^ 

However, unless £ becomes larger than l m no actual 
stripe correlations emerge. This happens only if, at 
even lower temperatures, the charge correlations are suffi- 
ciently strong to form a stripe liquid and different stripes 
are strongly correlated. We will focus on this temper- 
ature regime. Here, > r(T) > —2y/Q and thus G(x) 
exhibits an oscillatory behavior with £ > l m . 

It will be useful to introduce the positive dimcnsionlcss 
parameter e via 



4Q 



1 



(9) 



We then find (defining q 2 m = -r/2 ~ Q 1/2 with l m = 
that one can approximate the correlation function as 



S(q) 



(i) 



(10) 



In the last step we took only the leading term close to 
the peak at q m into account. By doing this approxi- 
mation we are breaking the charge neutrality condition 
(G{q = 0) = 0), only by a factor Q(q = 0)/G(q m ) = e 2 . 
For e <C 1 , the leading contribution of the Fourier trans- 
form of Q is given by 



G(x) 



c /£ sin(27nr/Z m ) 
Aire \x\ 



(11) 



which clearly shows the physical nature of the two length 
scales, £ = -=-, and l m = — . In order to have well corre- 
lated stripes, the correlation length £ must be much big- 
ger than the modulation length l m . This indeed trans- 
lates into e<1 which can be used as a small parameter 
of the theory. 

Within the large-A^ approximation, the temperature 
dependence of r is determined by 



r = r + u Q T 



8^ 



S(p). 



(12) 



For the case without frustration, Q = 0, the usual critical 
temperature = BZL_jp I results from the requirement 
r (T~) — 0. However, for finite Q no finite transition 
temperature occurs within the large- N approximation. 
Instead, one finds from Eq, 10 that 



{T)=ra+ 2^\2- +k 



(13) 



In Eq.[12|we ignored additional 1/N corrections, charac- 
terized by the self energy at zero momentum. In the limit 

o 2 

£ > l m under consideration we find T,q (0) ~ — ImS.^ 
Since e is small, the one loop self energy correction to G 
can be safely neglected, given that Eg is an additive cor- 
rection to the second term on the r.h.s. of Eqjl^ which 
behaves, in leading order, as e _1 . From Eq.| and Eqjl| 
we find at low temperatures 



2 W 



T/T c ° 



)l/4 



1 



T 



T ■ 



1 



(14) 



where in the last step |ro| ^> 2^/Q was assumed. This 
relationship will be useful for the determination of the 
temperature, Ta, where glassiness sets in. 



IV. THE STRIPE GLASS 

A. spontaneous ergodicity breaking and replica 
formalism 



With few exceptionscS~H3, the analytical investigation 
of glassiness due to the emergence of a large number 
of metastable states has concentrated on systems with 
quenched randomness. A major step forward was made 
in Ref.ElEU, where a new replica approach, equally appli- 
cable to quenched random and nonrandom systems, was 
developed. Within this approach, the configurational en- 
tropy for models of structural glasses was calculated, 
in good agreement with computer simulations £3 Here, 
we use this approach to investigate the physics of self- 
generated stripe glasses. 

The equilibrium's free energy density is given as F = 
— y logZ. It is of relevance only if the system is kinet- 
ically able to explore the entire phase space. Alterna- 
tively, one can introduce the averaged typical free energy 
of a system us ing t he following recipe of an "ergodicity 
breaking field"BEI . 

Locally stable field configurations can be identified us- 
ing a test field ip (r) and computing the partition sum 



Z[tp}= / Dtpe- n[v]/T -i 1 rf M</>(x)- v (x)] 2 



where g > denotes the strength of the coupling, 
dently, the free energy 



(15) 



Evi- 
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fbP]=-T log Z [ip] 



(16) with replicated partition function given by 



will be small when the field ip equals to a field configura- 
tion which locally minimizes TL. Thus, sampling all con- 
figurations of the V'-field, weighted with exp (— @f [tp]\ , 

is essentially a procedure to scan all locally stable con- 
figurations. The quantity 



F = lim ^ f DiP f m exp (~pNf ty] 



(17) 



is the weighted average of the free energy density of 
all locally stable configurations. Here, W = J Dip 

exp f [ip] /T^j is introduced for proper normalization. 

It is physically appealing then to introduce the free 
energy difference, 5F, via 



F = F-SF, 



(18) 



where 5F gives the amount of energy lost if the system 
is trapped into locally stable states and hence not able 
to explore the entire phase space of the ideal thermo- 
dynamic equilibrium. If the limit g — > on Eq. |l7j 
behaves perturbatively, 6F = 0. This indicates that the 
number of locally stable configurations stays finite in the 
thermodynamic limit, or at least grows less rapid than 
exponential with V. In this case all states are kineti- 
cally accessible. On the other hand, if the limit g — > 
does not behave perturbatively, it means that the num- 
ber of locally stable states, Af ms , is exponentially large 
in V. This allows us to identify the difference between 
the equilibrium and typical free energy as an entropy: 



SF = TS r 



(19) 



The configurational entropy density, S c — logA/" ms , is a 
measure of the number of metastable states and is an ex- 
tensive quantity if there are exponentially many of those 
states. Its emergence renders the system incapable of 
exploring the entire phase space. S c is then the amount 
of entropy which the system that freezes it into a glassy 
state loses due to its nonequilibrium- dynamics. 

In order to find an explicit expression for S c one intro- 
duces a replicated free energy:c3 



F (to) = — lim — log / D ipZ hp] 



from which F can be obtained as F = dm F( m ) 

am 

hence: 



S,= 



1 dF{m) 



T dm 



(20) 



and 



(21) 



Inserting Z [i/j] of Eq. into Eq. ^ and integrating over 
ip, one gets 



T 

F(m) = logZ(m) 

TO 



(22) 



Z(m) = lim J D m (pexpl~J2H [<p a ] /T 



(23) 



which has a structure similar to a conventional equilib- 
rium partition function. The ergodicity breaking field ip 
causes a coupling between replicas which might sponta- 
neously lead to order in replica space even as g — * 0. 
This order is then associated with a finite S c and thus 
glassiness. 

Formally, Eq.|23| equals the partition function of system 
with quenched random field analyzed using the conven- 
tional replica approach. The main difference is that, here, 
the limit m — * 1 has to be taken. The resulting many 
body problem in replica space is characterized by the 
matrix correlation function, Gab(l) = (^(l)^^ 
in replica space with Dyson equation: 



^" 1 (q)| 6 = ^ 1 (q)*-6 + s«6(q) 



g_ 

TO 



(24) 



Here, Go (q) is the Hartree propagator of Eq.|| which we 
approximate at low temperatures by Eq.|l0|. S a 6 (q) is 
the self energy in replica space. If we find that as a con- 
sequence of the ergodicity breaking coupling constant, g, 
S a b (q) has finite off diagonal elements we can conclude 
that there must be an energy landscape sensitive to the 
infinitesimal perturbation, g, supporting a glassy dynam- 
ics. On the other hand, if E a h (q) is diagonal, conven- 
tional ergodic dynamics results and the system is in its 
liquid state or may build an ordered solid. As pointed out 
above, this strategy is similar to the investigation of sym- 
metry breaking in conventional phase transitions where 
the off diagonal elements of an appropriately defined ma- 
trix self energy are associated with the order parameter 
of the transition (superconducting gap function or stag- 
gered magnetization in case of a superconductor or an- 
tiferromagnet, respectively). However, it turns out that 
in the present case the off diagonal elements of S a {, (q) 
jump discontinuously from zero to a finite value and a 
linearized theory with E a h (q) (1 — S a b) small (which de- 
termines the transition temperature in case of continuous 
phase transitions) will only give the trivial solution with 
vanishing off diagonal elements. A nonlinear theory for 
E a b (q) needs to be developed. 

Since the attractive potential between different repli- 
cas is symmetric with respect to the replica index we use 
the following Ansatz for the Green's function 



Gab (q) = (e(q)-^(q))<U+^(q) 



(25) 



i.e. with equal diagonal elements, Q (q), and equal 
off-diagonal elements, .F(q). Note, if one applies the 



G 



present replica formalism to systems with quenched dis- 
order it turns out that the replica symmetric Ansatz, 
Eq.^B], is equivalent to one-step replica symmetry break- 
ing in the conventional replica formalisms!]. The phys- 
ical interpretation of Q (r — r') = T _1 (ip(r)(p(r')) as 
thermodynamic (instantaneous) correlation function is 
straightforward. On the other hand, J- (r — r') = 
T _1 lim^oo (<p(r,t)tp(r' , 0)) can be interpreted as mea- 
suring long time correlations. As shown in the appendix 
A, inserting the Ansatz, Eq.^5], into Eq.^4| gives in the 
limit m — > 1: 



G- 1 (q) - Go 1 (q) + S e (q) 
for the diagonal elements, and 

-G 2 (q) (q) 



(26) 



^(q) 



1 - Q (q) E^ (q) 



1 



G~ X (q) - Ejf (q) 
= G (q) - K (q) 



(27) 



for the off diagonal elements, respectively. Here, Eg and 
Y^jr are the diagonal and off diagonal elements of the self 
energy in replica space. In the last equation IC denotes 
the deviations of the long time and instantaneous correla- 
tions. Analyzing the corresponding dynamical equations 
of the problem, it-Urrns out that IC is the static retarded 
response functionEa. In the liquid state fluctuation dissi- 
pation theorem gives IC — G and no long time correlations 
occur. 



The key assumption of the analytical approach to the 
self consistent screening approximation is that the off di- 
agonal self energy Ejr (q) is weakly momentum depen- 
dent. Specifically, we assume that due to the strong 
momentum dependence of the correlation function, the 
product (J(q)Ejr(q) varies with q predominantly due 
to £7(q). This assumption will be justified a posteriori. 
Also, our numerical results, which were obtained without 
any restriction on the q-dependence of Ejr(q), clearly 
show that this assumption is justified. We will then cal- 
culate E^r(g m ) at the modulation wave vector q m . 

>^From the analysis of the liquid state we know that 
G(q) is strongly peaked at the modulation vector q m with 
width By inspection of the Dyson equation, Eq,27, 
for T , one can see that for q ~ q m and E^-(g m )5(g m ) ^ 1 
(since G is peaked around q — q m ): 



Hi) < G{q). 



(29) 



Moreover, G vanishes rapidly away from the peak (as does 
Ei?(g m )G(g m )) and it follows from the same equation, |27], 
that for large \q — q m \ 



T{q) ~ -^{q)G 2 {q). 



(30) 



B. Defect wandering in stripe glasses 

The self energy-^n replica space was numerically in- 
vestigated in Ref.tS within the self consistent screening 
approximation which we summarize in appendix A. It 
was shown that below a characteristic temperature, Ta, 
an off diagonal self energy in replica space emerges, lead- 
ing to finite long time correlations, ^(q), as well as a 
finite configurational entropy density 

s c = S c /V (28) 

in the thermodynamic limit, V — > oo. Here we present 
an approximate but analytical solution of the same set 
of equations which has the appeal that the underlying 
physics of the stripe glass formation becomes much more 
transparent. It also reveals more directly the emergence 
of a new length scale which characterizes the wandering 
of defects in the stripe pattern after long times. In this 
context, we demonstrate that the melting of the glass as 
T becomes larger than Ta is a consequence of the fact 
that the characteristic length for defect wandering be- 
comes too large; the glass becomes too fluid causing the 
devitrification into a stripe liquid. 



If a solution for T exists, it is going to be a peaked at 
q m , but smaller and narrower than G- Consequently, if a 
stripe glass occurs, the long time limit of the correlation 
function in not just a slightly rescaled version of the in- 
stantaneous correlation function, but it is multiplied by a 
q dependent function that leads to a qualitatively differ- 
ent behavior for different momenta. Once a glassy state 
is formed, configurations which contribute to the peaks of 
G{q) and T{q), i.e. almost perfect stripe configurations, 
are almost unchanged even after long times. Close to 
q m T{q) is solely reduced by some momentum indepen- 
dent Debye- Waller factor exp (— D) = T{q)/G{q)- On the 
other hand, certain configurations which form the tails of 
G(q), i.e. defects and imperfections of the stripe pattern, 
disappear after a long time since now ^(q) <C G{q)- The 
ratio of both functions is now strongly momentum depen- 
dent. J-(q) becomes sharper than G(q) because certain 
defects got healed in time. Evidently, there must be a 
momentum scale (or equivalently a length scale) which 
determines the transition between these two regimes. In 
what follows we will identify and determine this length 
scale. A pictorial description of the defects on the stripe 
configuration and the meaning of A is presented on FigJJ. 



7 




FIG. 4. Pictorial description of the wandering of defects in 
the stripe pattern. The upper panel shows defects that can 
be healed by the wandering process. The lower pannel shows 
defects which are too far apart and cannot be healed. 

Due to our assumption that Ejf is weakly dependent 
on q, we concentrate on SjF (q m ) &t the modulation wave 
vector. One easily finds that E^r {q m ) < 0. A dimen- 
sional analysis furthermore shows that Hp is length -2 . 
This suggests to define a new length scale, A, via 



(31) 



For the subsequent calculation it is convenient to intro- 
duce in addition to the dimcnsionlcss parameter e which 
gives £ _1 = ^rp- a new dimensionless parameter, k, de- 
fined via 



A = 



Vk 2 - e 2 q n 



(32) 



Obviously, in the liquid state, where Ejr — > 0, we find 
A — > oo and it holds k = e. In a glassy state k > e. This 
Ansatz for Ejf, inserted into Eq.^7], yields 



K{q) = 



-1 

1- 



(33) 



for the correlation function K = Q — T . Note that K. 
has the same structure as Q but with e — > n. It imme- 
diately follows that it is the length scale A which deter- 
mines whether long time correlations are similar or dif- 
ferent from instantaneous ones. If \q — q m \ < A" 1 Eq.|29| 
holds, whereas for \q — q m \ > A -1 l ong time correlations 
are strongly reduced leading to Eq.|3(j. Consequently we 
identify A as the length scale over which imperfections of 
the stripe pattern manage to wander, i.e. defects can be 
healed, even in the frozen glass state. 

The next step is to determine Ejr(g m ) for a given value 
of A and to self consistently determine this length scale. 
The details of the calculation of E^r(g m ) with Q (q) as 
given in Eq.|l^ and K. (q) of Eq|3^ are summarized in the 
appendix B. The result is 



Ejr (q m ) 



7T 1 



1 



1 1 



(34) 



This has to be compared with our ansatz, Eq.pl|. To- 
gether with Eq.|3^ this gives 



*V (q m ) = - {k 2 - e 2 ) qt 



(35) 



Comparing Eqsj34| and |3^ we immediately find the non- 
linear algebraic equation 



7T 1 



1 



-Y 



(36) 



which determines the length scale A (via n) , as function 
of the correlation length and the modulation length, i.e. 
properties of the liquid state. One solution of this equa- 
tion is always e = K, which corresponds to the liquid 
state. Factorizing this trivial solution from [36] we find 
that the other solutions are given by 



TVS 



1 



1 



(37) 



This is a cubic equation which can be solved exactly. Be- 
fore we discuss this equation in some detail we analyze 
the condition for obtaining a nontrivial solution which 
corresponds to the onset of glassiness. For ^ ^> 1 the 
left hand side of Eqj3^ has its maximum ^ at n ~ 3e 
which gives the condition for the existence of a solution 
as — > 2. This defines the critical value of e as 



e A = e (T A ) 



1 

27T* 



(38) 



Thus, if the ratio of the correlation length and the mod- 
ulation length exceeds the critical value 



I n 



(39) 



T=T A 



long time glassy correlations emerge. As expected, be- 
sides strong frustrations, glassiness also requires suffi- 
ciently strong liquid correlations. Since it follows from 
Eq.|37| that n ~ 3e it also follows that 



A_ 

1 in 



T=T A 



(40) 



This type of behavior is evident on Fig.|| where, indepen- 
dent of the value of Q all curves for X/l m reach the same 
maximum value at T A . If defects and imperfections of 
the stripe pattern within the glassy state manage to flow 
over length scales larger than ~ |/ TO the glass becomes 
unstable because it is too fluid to support a frozen state. 
Thus, we identify the length scale A as the Lindemann 
length of the glass. 
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As can be seen on Fig.|], A is a monotonically increas- 
ing function of temperature. At small temperatures A 
grows linearly, evolving to a cusp at the dynamical freez- 
ing temperature, Ta- Above Ta, A becomes infinite and 
the devitrification is complete. 

Using the Q-dependence of e in the liquid state (see 
Eq.|l4|) the stripe glass - stripe liquid transition temper- 
ature is then given by 



SA 



J_ - yi/4 JWr c ° 

2n~2^ 1 T < 



which gives 



T A = 



7T 2 QV4 + 1 



(41) 



(42) 



where Tj? is critical temperature of the Q = prob- 
lem. Moreover, since the difference between Ta and Tk 
is small compared to (see inset of FigJ|), the Q de- 
pendence of Tk will be roughly the same as of Ta- As 
pointed out above, we do not expect Ta to be sensitively 
affected by an additional |<p 3 -term in the Hamiltonian. 
This is less clear for the Kautzmann temperature and the 
difference between Ta and Tk might well depend on the 
coupling constant v. Still we expect Tk to decrease for 
increasing Q. 




1.02 



FIG. 5. The Lindeman length A as a function of tempera- 
ture for different values of Q. Inset: Q dependence of Ta — Tk 



C. On lattice corrections 

Until this point our theory has been performed in the 
continuum's limit and effects due to lattice corrections 
have been neglected. The particular form of the prop- 
agator, Eqs.|| and [It], however gives rise to a specific 
sensibility of our results to lattice corrections which is 
worth mentioning. Within the continuum's limit, the 
large- A approach used in this paper yields no ordinary 



phase transition to a stripe-solid state, a result which 
has been pointed out earlierESEJ. This can be most eas- 
ily seen from the mean field equation, Eq.[l3|. A solution 
e = is only allowed if T = 0. For T > one always 
finds a large, but finite correlation length, £ = . The 
absence of an ordered state is a consequence of the dis- 
persion relation ui q , with 



A 



(43) 



In contrast to an antiferromagnet or a conventional 
charge density wave, where the low energy modes are de- 
termined by isolated points in momentum space, Eq.|43] 
gives rise to a d — 1 dimensional sphere of low energy 
modes in momentum space (see upper panel of Fig.^J) . It 
is this large phase space of low energy excitations which 
destroys long range order. 



,co(q) 





FIG. 6. Low energy modes. Upper panel: in the contin- 
uum's limit, A = 0, the phase space for low energy excita- 
tions is a d-1 sphere. Lower panel: The anisotropy gives rise 
to a "direction dependent mass" and the phase space for low 
energy excitations is reduced to a set of arcs 

The most dramatic effect of corrections beyond the 
continuum's limit is the appearance of an anisotropy on 
the low energy states. It is reasonable to assume that 
the Hamiltonian |l| has a next order correction of the 
type 4Kq 2 q 2 which leads to a direction dependent "mass" 
term in EqJ43] 



e 2 = e 2 {K = 0) + Kql 



sin 



(20 



(44) 



with angle <f> (in the x — y plane) and a dimensionless 
anisotropy parameter K . Of course this is only justified 
if K <§; q~ 2 . If however K ~ <7~ 2 the physics strongly 
depends on phenomena on the scale of the interatomic 
spacing. 

The mean field equation for finite K is then given by 



r(T) 



u T 



e (K = 0) 2 + Kql 



(45) 
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The system now undergoes a phase transition for arbi- 
trary small K at T c = 7tK - 1 / 2 a+i " ^ s ex P ec ted, T c van- 
ishes as AT — > . Moreover, as pointed out on Refs. 
T c does not merge with TJ? as Q — > 0. The origin of 
this phase transition is that the d — 1 dimensional sphere 
of low energy degrees of freedom is reduced to arcs of 
size dip ~ E/(K 1 ' 2 q m ) (see lower panel of Fig||). As e 
decreases, these arcs become indistinguishable from iso- 
lated points, a behavior similar to the case of an anti- 
ferromagnet or a charge density wave occurs, leading to 
an ordinary phase transition. The mean field analysis 
of the lattice version of Eq.?? of Ref.EI, which finds T c 
considerably smaller than T° supports K <C q^ 2 . 

In analogy to the mean field analysis of the liquid state 
one can also perform the theory of the glassy state for fi- 
nite K. If lattice corrections are strong and K ~ q~ 2 
a transition to a stripe solid occurs at T c . The low en- 
ergy excitations are located at isolated points and the 
behavior is equivalent to the one of an unfrustrated sys- 
tem. In this case, the glassy state will only occur if the 
solidification is avoided by supercooling. On the other 
hand, if K -C g~ 2 the low energy modes are unchanged 
and the lattice corrections become irrelevant. Since the 
glass transition does not require e to vanish, but solely 
to reach a certain finite limit ~ J- , we conclude that for 

„-2 

K < the glass transition is essentially unchanged. 



D. The configurational entropy 

The main argument for the emergence of a glassy 
state is the occurrence of an exponentially large num- 
ber if metastable states, characterized by the configura- 
tional entropy. S c is determined from Eq.^l] from the 
replicated theory defined in Eq.p3|, which gives F (m) — 
— — log Z(m). Within the self consistent screening ap- 
proximation it follows 



F (m) / (2mT) = tr log Q~ 



trlogP" 1 - trSg 



Since all quantities are matrices in replica space with a 
structure given in Eq.^5] the evaluation of expressions like 
trlogCf -1 etc. becomes straightforward. Performing the 
derivative with respect to the number of replicas accord- 
ing to Eq.Ell gives immediately for the entropy density: 



,(2) 



with the two contributions: 



and 

«(2) : 



2 J { g(q)J 



2 J 8vr 3 \ n \ l + « n e (q) 



<?(q) 



v Tljr(q) 



For the definition of II^-, Hg etc. see appendix A. Using 
the same approximations as for the evaluation of the self 
energy, Ejf, in appendix B we find 

^(q) k 2 - e 2 



0(q) 



The evaluation of the integrals is straightforward and we 
find 



,(1) _ <lrnK 
47T 2 



(46) 



Obviously, s c ^ only if k > e i.e. the Lindemann 
length, A, is finite. s c vanishes in the liquid state where 
k = e. The results for s c (T) for different values of Q are 
given on Fig.^ 

Using the results of the previous section where we 
found that for T = Ta the dimensionless quantities e 
and k take fixed values, we obtain at T = Ta- 



S c (T A ) = CVQ 



3/4 



(47) 



withC= (l + J^-L + ^md) = 1.1816X10- 3 . 
The configurational entropy decreases for decreasing Q. 
Since the modulation length behaves as l m ~ Q 1 ^ it 
follows S c oc The larger the modulation length, 

the smaller is the number of states one can form which 
clearly demonstrates that locally stripe correlations stay 
intact in all these configurations. It is more the packing 
and arrangement of defects which distinguishes different 
metastable states. Less of those packings are possible per 
unit volume if the modulation length grows, for simple 
geometrical reasons. The relation, Eq.^ was recently de=. 
rived by us using the concept of replica bound states.EZl 
The fact that we obtain the same result using an entirely 
different approach to solve the problem increases our con- 
fidence in the applicability of the self consistent screen- 
ing approximation, which also allows us to calculate the 
constant C. We can also compare S c (Ta) with the, nu- 
merical results of Ref.E3. S c was calculated in Ref.ll3 for 
the values Q — 0.01 and Q — 0.001, where we found 
S c = 3.45 x 10~ 5 and 6.4 x 10~ 6 , respectively. ^From 
Eq.(47| we find the values 3. 7 x 10~ 5 and 6. 6 x 10~ 6 , 
which agrees well with the numerical results. 

Finally, the condition S c (T = Tk) = determines the 
Kauzmann temperature below which no entropic advan- 
tage of the liquid state, compared to the amorphous solid 
state (the glass) exists. Even if one manages to anneal 
the liquid down to Tk without freezing into a glass, some- 
thing which might me achieved using an infinitely slow 
cooling rate, a mandatory transition into an amorphous 
solid occurs at Tk- In Fig.|^ we show the temperature 
dependence of S c for different values of the frustration 
parameter Q as well as the Q dependence of the slope 
(dS c /dT) TK . 
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T/T k (Q) 



FIG. 7. The configurational entropy for different values of 
Q. Inset: dependence of (dS c /dT) T on Q 



E. Dynamics and flow via entropic droplet formation 

So far, we have shown that within the self consistent 
screening approximation of Eq.jl], an exponentially large 
number of metastable states occurs. The conclusion that 
somewhere below Ta nonequilibrium dynamic sets in is 
actually obtained using a purely thermodynamic char- 
acterization of the spectrum of metastable states. It 
is very plausible that an exponentially large number of 
metastable states is necessarily connected to glassy dy- 
namics. In several random spin models as well as in the 
model of self generated glassiness in frustrated Josephson 
junction arrays, this point of view has been clearly sup- 
ported by actual dynamical calculations. For so called 
infinite range models (i.e. within mean field approxima- 
tion), where the barriers between the various roatastable 
states diverges, freezing at Ta has been found.OC3 

Solving the Langevin equation for the time evolution 
of the correlation and response function of cp (x, t) , us- 
ing the supersymmetric formulation of the Martin-Siggia- 
Rose approachE3 within the self consistent screening ap- 
proximation, we find that the emergence of exponentially 
many metastable states also .leads to stripe glass state 
within a dynamical approach£3 Interestingly, this mode 
coupling type approach gives exactly the same criterion 
for the emergence of glassiness as the above replica ap- 
proach only if one properly takes the aging behavior|-af 
the dynamical evolution into account, following Ref.c3. 
Thus, the replica approach employed here takes the ef- 
fects of aging correctly into account. 

Within the mode coupling or replica integral equation 
approaches a perfect freezing occurs at the temperature 
Ta- However, more realistically, Ta is rather a crossover 
scale to a regime with slow activated dynamics and not 
the actual freezing temperature. Depending on the his- 
tory of the system, the laboratory glass transition where 
the time scale of motions exceeds a certain limit occurs 
somewhere between Tk and Ta- A key question in 
this context is the nature of the dynamical processes for 



T K <T< T A - 

This dynamics involves droplets or instantons, essen- 
tial singularities from the point of view of perturbative 
approaches. A complete formal treatment is therefore 
difficult. Nevertheless, a reasonable description of the 
dynamics in a system with finite S c is given in Ref.O, 
where the nucleation of droplets with size R of a new state 
within an old one was argued to be the main dynamical 
processes. The free-energy gain of a droplet formation 
is caused by the entropic gain due to the exploration of 
new states, i.e. by Ts c i? 3 , whereas one has to take into 
account that such a droplet implies a finite surface en- 
ergy, characterized by a scale dependent surface tension, 
a (RX, A renormalization group calculation, based on 
Ref.cj, leads to the size dependent surface tension 



a(R) = cr (RAy 



(48) 



with 9 — reflecting the fact that the interface be- 
tween two states is wetted by intermediate states. This 
analysis leads to an characteristic energy barrier AE oc 
(Ts c (T))~ 1 which implies a characteristic relaxation time 
which follows a Vogel-Fulcher law 



r oc exp 



DTjc 
T - T K 



(49) 



independent on the dimension. Here, fragility parameter 
of the Vogel-Fulcher law is given by 



D 



3^o 2 



T 2 T, 



K 



dS c 
dT 



(50) 



Xia and Wolynesa have shown that this scenario gives a 
quantitative description of viscous flow m molecular liq- 
uids A straightforward extension of Refoo along the 
lines of Ref.cZl also gives a width of the distributions of 
different activation energies characterized by the mean 
square width (<5i? 2 ) of the droplet size, which might be 
compared with the results of Ref.EJ. This droplet picture 
implies that the glass breaks up into domains of differ- 
ent metastable states, separated by wetted surfaces (con- 
sisting of intermediate states) leading to a rather small 
surface tension. This physical picture is very similar to 
what is usually called a " clusteE-,spin glass", motivated 
by the observations made in Reffl based on NMR exper- 
iments. 

For a quantitative analysis we have to develop a the- 
ory for the bare surface tension er for the stripe glass. 
This differs from the theory of molecular liquids because 
of the long range forces in the model. In what follows we 
give simple estimates for ao , based on a variational argu- 
ment. There are two sources of the bare surface tension 
in Eq.|l|, the gradient term and the long range Coulomb 
term. Assuming a droplet configuration with locally or- 
dered charge configuration 



(p (r) = (po cos (2itr /l m ) tanh 



R 



(51) 
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with droplet radius, R, and wall thickness, l w . In order 
to make progress, we assume that the thin wall limit, 
l w <C R applies and find for the gradient term: 

EP = inipll- 2 [ R+ W ' r 2 dr ~ ^l^R 2 . (52) 



and for the long range Coulomb term 

~ 4ny 2 l- l R 2 . (53) 



For larger values of the frustration parameter, Q, l w is 
estimated by the largest length scale of the problem (ex- 
cept R of course) which should correspond to the lowest 
energy of the droplet wall. This gives l w = £. Alter- 
natively, for smaller values of Q the surface tension is 
dominated by the contribution of the Coulomb term and 
we find 

CT = 47T^/~ 1 . 



^From the same variational argument we find tp 2 = 
-^Q 1 ^ 2 , which gives, together with l m ~ 2nQ 1 l A , the re- 
sult (To — §Q 3 ^ 4 - Thus, if the frustration parameter de- 
creases, the barrier height between different metastable 
states disappears. Even though the number of metastable 
states an thus ^ \ T decreases for decreasing Q, the sur- 
face tension term dominates and it follows that 



D (Q -> 0) -> 0. (54) 



The fragility parameter does not vanish according to 
a power law, mostly because Ta — Tk, which enters 
It — ta-tL ^ as logarithmic behavior as Q — > 0. 
Essentially, 13 vanishes with Q similar to a square root. 
In Fig.|| we compare D (Q) with the results obtained in 
Monte Carlo sipulations of the lattice version of Eq.[l] by 
Grousson et a/.c2l. Here we multiplied the result from Eq. 
|50| by an overall prefactor, leading to a very good agree- 
ment between our analytical theory and the numerical re- 
sults oc3. Since the calculations in Ref.E3 are performed 
for a lattice version of Eq.|l| the actual T dependence of 
the correlation length differs from ours. Correspondingly 
the absolute magnitudes of Tk, which enters in Eq. |5^, 
are different. For this reason we do not expect the ab- 
solute magnitude of D, but solely its Q dependence to 
agree in both approaches. 
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FIG. 8. Comparison of D with the results obtained in 
Monte Carlo simulations of the lattice version of Eq.[l] by 
Grousson et al. 



V. CONCLUSIONS 

Glasses are typical examples of systems of many in- 
teracting particles that have a tendency to self-organize 
into mesoscopic structures. In this paper we studied 
the slow activated dynamics of charge inhomogeneities in 
doped Mott insulators. We developed an analytical ap- 
proach which enabled us to identify the underlying phys- 
ical mechanism for glassiness in a uniformly frustrated 
system. We showed that when the charge correlations 
are sufficiently strong, specifically, if f- > 2, the stripe 
liquid - stripe solid transition can become kinematically 
inaccessible because the system undergoes a glass tran- 
sition, driven by the emergence of an extensive configu- 
rational entropy. We demonstrated that at this point a 
Lindemann length, A, emerges, which is a length scale 
over which imperfections of the stripe pattern manage to 
wander. Finally, we apply our results to the scenario of 
Ref.cJ to calculate the characteristic relaxation time of 
the nonequilibrium state. We concluded that the charge 
fluctuations in doped Mott insulators have a tendency 
to self organize into droplets of metastable states, dis- 
tinguished by the packing of stripes with different orien- 
tation and the arrangement of defects. These droplets 
relax according to a Vogel-Fulcher law, characteristic of 
structural glasses. We further compare our results with 
the doping dependence of Tk, as deduced from experi- 
ment in Refo, and show that it is properly described 
within our theory. 
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APPENDIX A: SELF CONSISTENT SCREENING 
APPROXIMATION 



In this appendix we summarize the self consistent 
screening approximation which was the basis, of the nu- 
merical investigation of stripe glasses in Ref .113 and which 
is the framework in which we determined the diagonal 
and off diagonal element in replica space of the self en- 
ergy, leading in particular to Eq.[}4| The details of the 
calculation of Eq.[?4| are then given in appendix B. 

Eq.^3] has a formal similarity to the action of the ran- 
dom field Ising model, obtained within the conventional 
replica approach, whicb-allows us to use techniques, de- 
veloped for this modelE-X Introducing an A-component 
version of Eq.[j] with field ip — (ipi, (Pn) and coupling 



constant, u 



- ft, with 



ted «o we use a self consistent 
screening approximatioro, which is exact up to order 
1/A. At the end we perform the limit N — 1. The appli- 
cability of this approximation is supported by the strong 
indications for a-»on-singular large-A limit of Eq.|l|, as 
discussed in RefBI. 

Before we discuss the self consistent screening approx- 
imation, it is useful to summarize a few properties of 
matrices in replica space with structure similar to |2^. 
Introducing a matrix E such that E a fc = 1 and the unit 
matrix 1, it is easy to see that the product of any two 
to x to matrices with structure 



A = ail + a 2 E, 



(Al) 



is given by 



AB = (ai£>i) 1 + (0462 + ^201 + toq^^) E. 



This leads to 



«2 



ai ai (ai + ma 2 ) 



-E 



(A2) 



for the inverse of A. This property was used for the 
derivation of Eqs. Gfl and and will be used below. 



+ 



FIG. 9. Diagrams for the self consistent screening approx- 
imation 

The self consistent screening approximation is de- 
scribed by the set of Feynman diagrams shown in FigJ|. 
The self energy is given as 



Sab (q) 
where 



v / V ab (p)G ab (p + ci). (A3) 



Z>(P) 



n(p))" 



(A4) 



is determined self consistently by the polarization func- 
tion 



Hah (p) = / Gab (q + P) Gba (q) 

J (2tt) 



(A5) 



In the above set of equations the p-integration has to 
be cut-off at |p| = A and the temperature, T, and the 
coupling constant, uq, occur only in the combination 
t>o = uqT. The Ansatz ^5] for the Green's function im- 
plies an analogous structure for E a fc (q) and H a b (q) in 
replica space. Inserting this ansatz into n ao (p) gives 



n = (n e -LV)i + n^E 



(A6) 



where the diagonal and off-diagonal elements of the po- 
larization function are 

d 3 q 



n e ( P ) = J ^3S(q + p)S(q) 



^(P) = J ^3^(q + p)^(q). 



(A7) 



Using the rule |A2j, it is now straightforward to deter- 
mine T> a b (p) which leads, in the limit to — > 1, to 



where 



V = (Vg - Vjr) l+Vjr'E 

Vg (p) = {vo 1 +U g (p))" 1 



(A8) 
(A9) 
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and 



?Mp) 



n^( P )2?§ ( P ) 
i-n^ ( P )v g ( P ) 



(A10) 



Analogously, inserting the above equations into A3 
get for the self energies: 



we 



S^) 1 + EjfE 



where 



S S (q) = | / 7^3^ (P) G (P + q) (All) 



and 



5V (q) 



(2tt) 

d 3 p 



NJ (2nf 

The set of equations is closed by the Dyson equation, 

E 

Q- 1 (q)L = (Go 1 (q) + £g - Sjf) i + s^e. (A12) 

that give s, according to in the limit m-> 1 the equa- 
tions]^ aud 

In Ref.ta it was shown that this coupled set of equa- 
tions gives Yjjr (q) ^ below a characteristic temperature 
which was related to the occurrence of glassiness. In the 
next appendix we present an approximate analytical so- 
lution of this problem. 



APPENDIX B: ANALYSIS OF THE SELF 
ENERGY 

In this appendix we give the main technical details for 
the calculation of the off diagonal self energy in replica 
space, T,jr. The diagonal self energy T,g, which turns out 
to be negligibly small compared to the leading mean field 
terms, can be determined in, a ..very similar fashion and 
was already analyzed in Ref.E£rc3 Sjf and Eg were cal- 
culated numerically within the self consistent screening 
approximation in Ref.Ej. The virtue of the analytical 
calculation presented here is that it is much more trans- 
parent. In all steps of this calculations we did check 
the reliability of our approximate analytical treatment 
by comparing it with the numerical results. 

We start by calculating the polarization function 
n 6 (g): 



u c{q) = J ^§£p£ P +q 



= < tan 



q 

eq„ 



1 

— tan 

2 



— tan 
2 



2<Zn 



sqn 



ggm - q 

eq m 



1 „ ^ 2eq„ 

87rg m e ;1 y ^ it 

B(2q m -q) 2eq„ 

TEq~p q ^ IT 



where we used the approximate expression, Eq.[10], for the 
correlation function Q (x). An analogous calculation for 
the off diagonal polarization function Il^-(g) gives 



|-§/C (p) K (p + q) 



as well as 



0s(p)K(p + q) 



1^ 2^ 

<S>{2q m -q) 2 K q„ 
16n 2 q q ^ TT 



q < ^qn 



87rq m K 3 
e(2q m ~q) K+£ „ 



For (e, K)q m < q < 2q m , we can use the approximate 
expressions 



n e (g) 



1 



16qe 2 



Ibq \e k 



which gives 



Dg(q) = 



I6q m a 

1 + 2^4 

q s z 



with dimcnsionlcss number a = , „"° < 1. a vanishes as 
0, but it always holds that e 2 <§; a. Note, for the 



T 



numerical solution in Ref. it holds a ~ 0.4. Combin- 
ing these results and using the fact that e 2 /a <C 1, the 
product Dg{q)Yljr(c[) becomes momentum independent: 



D g (q)ILjr(q) ~ 



and, as a result, Dyr(q) becomes proportional (by a factor 
smaller than 1 in magnitudes) to Dg(q) 



D F {q) = 



-Dg (g)ILjr (q) 



1 - Dg (q) (q) 



Dg (q) 



l-(l-f) 



Dg(q) 



We are now in the position to analyze the self energy 
S^: 



^(q) = / 0^(q + p)^(p) 



(Bl) 



Since D?r (q) is only weakly momentum dependent the 
same holds for Sjf (q) and we can estimate it at the mod- 
ulation vector. It follows with t = \J2 (1 + cos#) 



5V (<Zm) 



tdtD T (q m t) I ^FiP) 
p 2 dp , 



Dt (q m ) 



4tt : 



(B2) 
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Using T{p) = G(p) - £(p) and J (p) = 

_L(fe +A ) aswellas/^/Cb) = &(^ + A), we 
find 



E.F (g m ) = -- 



(1 



1 



(B3) 



Thus, we have determined the off diagonal self energy at 
the modulation wave vector as function of the three es- 
sential length scales of the problem: l m = £ = 

and A = . , 2 — . Note, the dependence on the mo- 

VK—sq m 

mentum cut off, A, cancels completely making this result 
robust against lattice corrections. 



APPENDIX C: COMPARISON WITH THE 
FRUSTRATION LIMITED DOMAIN SCENARIO 

The possibility of glass formation of the model, Eqjl], 
has been pointed out in Ref.Ea. As discussed in this ap- 
pendix, we disagree with the detailed argumentation of 
Ref.Bj. However the recognition that a model of the kind 
of Eqjl] can potentially describe glass formation was a 
very important observation. 

The aim of Ref.E£l was to present an alternative sce- 
nario for glassiness in structural glasses formed of under- 
cooled molecular liquids. Though the microscopic justi- 
fication of Eq.[l] for the description of structural glasses is 
at the least unresolved, one can yet, in principle, imagine 
that such long range interactions are caused by a scenario 
like that based on icosahedral order which is frustrated 
by the lack of Euclidean curvature of the effective space.B 
In what follows we will solely consider Eq.[j] as a given 
model and leave aside whether it applies to stripe glasses 
in doped Mott insulators las we claim) or to molecular 
liquids (as claimed in Re£,3). 

The main idea of Ref.ca is that due to the frustrating 
interaction of Eq|l| the system is broken up into ordered 
domains of size Rd ~ £,~ 1 Q~ 1 ^ 2 with £ being the correla- 
tion length that controls the fluctuations inside ordered 
domains. Furthermore, £ <C Rd is assumed. Within 
each domain the ordering essentially corresponds to the 
one of a finite Q = system (not of a system with stripes 
as in our entropic droplet picture). Since for Q = 0, £ 
diverges at the critical point , the avoided frustration 
scenario suggests that for Q > 0, £ (T) still has a peak 
close to T° only rounded due to the frustration. It is 
then argued that there is a relaxation rate r -1 due to 
the reorganization of domains which obeys an Arrhenius 
form 



r" 1 cx e -AF(fl D )/T 



(CI) 



where AF (Rd) ~ f is the activation energy of 

the domain. Finally, it is asserted that the divergence in 
the viscosity of the system at low temperatures is deter- 
mined by r\ cx r. 



We disagree with this picture. As shown below, the 
scale Rd, which signals the relevance of the long range 
interaction, can easily be identified as the Debye screen- 
ing length of charged particles of size £ with the expected 
charge density. Even when the Debye picture applies, 
we find it hard to understand how conventional screen- 
ing can lead to a dynamics which is dominated by the 
activated reorganization of screening clouds and where 
the natural Langevin description gives a fast relaxation. 
When £ <C Rd there are many short wavelength excita- 
tions of weakly coupled charges and the system should 
rather behave as a high temperature plasma than, as a 
glass. In addition, for N — > oo it was shown in Rcf.ES that 
Debye screening occurs only at temperatures T > T®. 
At low temperatures modulated structures result which 
lead to the formation of well correlated stripes as £ > l m . 
£ (T) does not exhibit a maximum at a temperature com- 
parable to Tj: but keeps growing until either a stripe solid 
or a stripe glass is formed. 

We now show that the scale Rd ~ £ _1 Q -1 ^ 2 is the 
Debye screening length of charged particles of size £. We 
perform a coarse graining of the system into regions of 
linear size £ centered around positions X^. Then Eq.|l| 
becomes: 



i>j 



(C2) 



where 



Hi ~ \ J a d 3 x {rw(x) 2 + (V^(x)) 2 + ^(x) 4 } (C3) 



and charges 




^- / d 3 xtp(x) oc VQC 



5/2 



In the last step we assumed that within the volume oc £ 3 
the system is essentially ordered and used tp oc {^~PI V 
with the critical exponents j3 = i and ^ = (d — 2) 4 
obtained within the large N ap proximation, for d = 3.E2I 
The usual analysis of Eq.C2 within the Debye approx- 
imation, i.e. solving the Poisson equation with induced 
charges distributed with Boltzmann weight, gives the De- 
bye screening length: 



,-2 _ 27rg 2 

l D " rp 



-n, 



where n is the density of charges, here given by n ~ £ 3 . 
This finally gives 

Id - Q- 1/2 C\ 

which is precisely the length Rd proposed on Ref.B 

In summary, it is unclear why Debye screening should 
cause a break up of the system into ordered domains of 
size Rd = Id (remember £ -C Id) an d a slow activated 
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dynamics of such domains, necessary to obtain large vis- 
cosities. Second, for the scenario of Ref.CJ to work, one 
also has to assume that the correlation length decreases 
at low temperature and that Debye theory still applies. 
The actual analysis of Eq.[l] for M — > oo does not show 
both assumptions to be justifieaEl However, it is inter- 
esting that for N < oo the emergence of a Josephson 
length scale at low temperaturescj might give rise to a 
competition of physics on two distinct length scales which 
then, in principle, could lead to a new long time relax- 
ation along, the lines of the frustration limited domain 
approach.til It is interesting to explore the relationship 
of this scenario to our replica based theory. 

Finally, it is worth pointing out that the entropic 
droplets discussed in this paper are qualitatively different 
from the domains introduced in the approach of Rcf.E3. 
Whereas the latter correspond to thermodynamically sta- 
ble configurations, similar to domains in ferromagnets 
caused by the long ranged dipole-dipole interaction, our 
entropic droplet are formed by the various metastable 
states. Transitions between different droplets are caused 
by a gain of entropy of a system in what would be an 
otherwise frozen noncquilibrium state. 



17 J. Schmalian H. Westfahl Jr. and P. G. Wolynes, to appear 
on Int. Journ. of Mod. Phys. 
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